Dynamics of Impurity and Valence Bands in Gai^Mn^As 
within the Dynamical Mean Field Approximation 
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We calculate the density-of-states and the spectral function of Gai-ajMn^As within the dynamical 
mean-field approximation. Our model includes the competing effects of the strong spin-orbit cou- 
pling on the J = 3/2 GaAs hole bands and the exchange interaction between the magnetic ions and 
the itinerant holes. We study the quasi-particle and impurity bands in the paramagnetic and ferro- 
magnetic phases for different values of impurity-hole coupling, J c , at the Mn doping of x = 0.05. By 
analyzing the anisotropic angular distribution of the impurity band carriers at T — 0, we conclude 
that the carrier polarization is optimal when the carriers move along the direction parallel to the 
average magnetization. 
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PACS numbers: 75.50.Pp,2.70.-c,5.10.-a,71.27.+a 



The combined magnetic and semiconducting charac- 
teristics of ferromagnetic semiconductors make them ex- 
cellent candidates for spintronic applications 0. In par- 
ticular, GaAs doped with Mn (Gai-^Mn^As) is promis- 
ing as a spin-carrier injector in spintronic devices 0] due 
to its relatively high magnetic transition temperature Q 
and its potential integration within the current semicon- 
ductor technology. However, the properties of magnetic 
semiconductors still need to be greatly optimized since 
practical uses demand magnetic transitions above room 
temperature and carrier polarization of at least 70% 0. 

A realistic model which incorporates the relevant 
bands of the host material is crucial to guide the exper- 
imental efforts in the search for optimal magnetic semi- 
conductors. In Gai-^Mn^As, the Mn ions are in the 
Mn 2+ state with a half-filled d shell of total spin S = 5/2 
0, 0. Since Mn 2+ ions primarily replace Ga 3+ , they 
contribute carrier holes to the p-like valence band. The 
strong spin-orbit interaction couples the I — 1 angular 
momentum to the electron spin (s = 1/2), resulting in 
a total spin J = \l + s| = 3/2 for the two upper va- 
lence bands and J = \l — s\ = 1/2 for the split-off band 
0. Since the J = 3/2 bands are degenerate at the F 
point, an accurate model should include at least these 
two bands. However, a more realistic approach should 
incorporate the split-off and conduction bands as we dis- 
cuss later. 

Here, we continue our dynamical mean-field approxi- 
mation (DMFA) HHHII3 study of the effects of strong 
spin-orbit coupling in Gai-xMn^As^lJ. While we previ- 
ously examined the influence of the spin-orbit interaction 
on the ferromagnetic transition temperature T c and the 
carrier polarization [Tlj , we now focus on the density-of- 
states, the spectral function, and the dispersion of the 
quasi-particle and impurity bands. We also discuss the 
anisotropy of the spectra in the ferromagnetic phase and 
its influence on the transport properties. 

Although the formation of the impurity band has been 



captured in previous DMFA studies ^2] , their model does 
not take into account the spin-orbit coupling and is un- 
able to address the reduced carrier polarization within 
the impurity band. The DMFA describes the impurity 
band through quantum self-energy corrections which are 
not included in other mean-field theories. Because this 
method is non-perturbative, it allows us to study both 
the metallic and impurity-band regimes as well as both 
small and large couplings. Although the precise role 
played by the impurity band in Gai-^Mn^As is still 
controversial, an array of experimental probes, such as 
ang le-resolved photoemission |13|| . infrared spectroscopy 
|l4lll5lll6| . spectrosco pic ellipsometry |17|. scanning tun- 
neling microscopy Il3 , and photoluminescence tech- 
niques po| , display features characteristic of an impurity 
band. 

We start with the Hamiltonian proposed in Ref . fill and 
ED: 



H — H — J c ^ Si • J(Ri 



(1) 



The first term incorporates the electronic dispersion and 
the spin-orbit coupling of the J = 3/2 valence holes 
within the spherical approximation . The second term 
represents the interaction between the Mn spins and the 
valence holes [23|, with J c the exchange coupling and 
J(i?i) the total J = 3/2 spin density of the holes at the 
site i of a Mn ion with spin S^. The relatively large mag- 
nitude of the Mn spin (S — 5/2) justifies its classical 
treatment. 

Since the typical hole concentration is small (around 
5%), the holes gather at the J=3/2 bands around the T 
point. This supports the use of the spherical approxima- 
tion |22j, for which the non-interacting Hamiltonian of 
pure GaAs is rotationally invariant. Hence, Hq is diag- 

— k 2 



onal in a chiral basis, Hq — 



^ 2m. 
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k 7 Ck 7 , where c k 
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creates a chiral hole with momentum k parallel to its 
spin and J • k = ±3/2 or ±1/2. The two band masses 
« 0.5m and m; ps 0.07m correspond to the heavy 
and light bands with 7 = ±3/2 and ±1/2 respectively 
(771 is the electron mass). Notice that for convenience 
we use the hole picture so that the valence bands have 
a minimum instead of a maximum at zero momentum. 
However, in displaying the results we reverse the sign 
back to accommodate the usual convention. 

As discussed previously ^lj, the coarse-grained Green 
function matrix in the non-chiral fermion basis is 



1 



G(iu n ) = }^[iuj n I - e(k) + fil - E(iw n )]- X , (2) 

k 

where N is the number of k points in the first Bril- 
louin zone, fi is the chemical potential and e(k) = 

. - k 2 - - 

PJ(k)- -R(k) is the dispersion in the spherical approx- 
imation. Here, R are spin 3/2 rotation matrices that re- 
late the fermion operator Ck 7 to its chiral counterpart 
Ck 7 = i? 71 ,(k)ckj/, and repeated spin indices are summed. 
The mean- field function G(iio n ) = [G (iu> n ) + Ti(iu/ n )] 
is required to solve the DMFA impurity problem. At a 
non-magnetic site, the local Green function equals the 
mean- field function G non = Q, while the local Green 
function at a magnetic site is G mg (iu> n ) — [Q (iuj n ) + 
J C S • J] . To obtain this result, we treat disorder in a 
fashion similar to the coherent potential approximation 
(CPA) [p. 

Now G mg (iw n ) must be averaged over all possible spin 
orientations at the local site and over all possible impu- 
rity configurations on the lattice. The former is imple- 
mented by introducing the angular distribution function 
exp[-5 e //(S)] 



P(S) 



where Z = J dfls sxp[S e ff(S)] 



and S e jf(S) is the effective action of the system [25j 



Seff{S) = -^logdet ^(iWnX^ -1 («■>») + J C S • J) 



(3) 



The extra factor of G(iu) n ) in Eq.J3J) is introduced to aid 
in convergence. If the Mn ions are randomly distributed 
with probability x, then the configurationally-averaged 



G mg (iLO„))x + 



Green function reads G avg (iuj n ) = 

G(iu n )(l - x). 

When the magnetic order is along the z axis, G mg , Q 
and S are diagonal matrices and the angular distribution 
function depends only on the polar angle of the impurity 
spin: P(S) = P{0) |26|. In fact, G mg can be written 
as a combination of the identity matrix, I, J z , J| and 
Jf , where the coefficients are of order zero, one, two and 
three, respectively, of the magnetization: G mg — 0(1)1+ 
0(M)J z + 0(M 2 )J| + 0(M 3 )jf 0. Making use of 



these results a new algorithm has been designed which 
dramatically reduces the computational time compared 
to our earlier work ^l]]. This improvement is due to 
the fact that the rotational symmetry of the model is 
explicitly broken by our choice of a preferential direction 
for the magnetization. Also, since the symmetry breaking 
is already incorporated in this algorithm, we no longer 
need a small magnetic field to break the symmetry along 
a magnetization axis [TH . 

To calculate dynamical quantities, we work in the real 
frequency domain, where the coarse-grained Green func- 
tion matrix is 



e(k)-E( W )]- 



(4) 



with Q = lu + i0 + . 

Generally the Matsubara and real frequency Green 
functions need to be iterated simultaneously and aver- 
aged over the polar angle distribution P(9). However, if 
we focus only on the paramagnetic phase (T > T c ) and 
the ferromagnetic ground state (T = 0) the angular dis- 
tributions are known and self-consistency in the Matsub- 
ara domain is not necessary. In the paramagnetic phase, 
the angular distribution of the Mn spins is completely 
random so that P(0) = 1/tt. In the T = ferromag- 
netic ground state, the average impurity magnetization 
achieves its full value and P(8) — 6(9). 

After the coarse-grained Green function is self- 
consistently calculated, the density-of-states is computed 



DOS(lo) 



1 



Im Tr G(Q) , 



(5) 



where Tr is the trace. Each of the diagonal elements 

of Im G(f2) is the projection of the density-of- 

states onto a state with fixed J z component, i.e. J z — 
+3/2, +1/2, -1/2, and -3/2. 

We are also interested in the spectral function defined 

as 



A(k,u) 



1 



Im Tr [QI - e(k) - t(uj) 



(6) 



The center of the quasi-particle peak in the spectral func- 
tion represents the renormalized quasi-particle energy 
w M (k) (f-t = 1,2,3,4), which can be obtained by solving 
the condition: 



Re [ni - e(k) - ±(uj)} diag = 



where the subscript "diag" means that we first diagonal- 
ize the matrix, then solve the equation for each diagonal 
element. 

We focus on the doping of x = 0.05, for which T c is near 
the highest reported 0, H El |H |3fl HJ |H H3 . Fig.ffl 

shows our results for the density-of-states for various val- 
ues of J c in the ferromagnetic phase (main panel) and the 
paramagnetic phase (inset). As J c increases, states with 
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FIG. 1: Density-of-states for J c =1.0, 2.0, 3.0, and 4.0 eV, at 
T — (main panel) and T > T c (inset). The appearance of 
an impurity band is evident for J c > 2.0 eV. The broader im- 
purity band at T = 0, compared to that at T > T c , indicates 
that lifetime of the bound-state holes decreases as the system 
becomes magnetically ordered. 
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FIG. 2: Decomposition of the density-of-states in terms of 
J z components, for J c — 3.0 eV at T — (main panel) and 
T > T c (inset). The chiral nature of the holes due to spin- 
orbit coupling mixes up states with different J z components, 
making the total DOS only partially polarized at T = 0. 



positive energy appear inside the semiconducting gap. 
These states correspond to the Zeeman splitting of the 
hole levels induced by the local impurity magnetization. 
For J c > 2.0 eV an impurity band clearly appears. For 
J c = 3.0 eV and T = 0, a second impurity band has 
started to form though it has not yet separated from the 
main band. We also observe a second impurity band ap- 
pearing at positive energies in the paramagnetic phase 
when J c > 5.0 eV (not shown in the graphs). The ap- 
pearance of two impurity bands is consistent with the 
fact that the model includes two bands with J ■ k = ±3/2 
and ±1/2. 

As expected, the center of the impurity band shifts 
to higher energy as J c increases. However, our predic- 
tions for the energy of the impurity band are too large. 
We believe that this is a consequence of excluding the 
conduction band from our model, since band repulsion 
with the conduction band pushes the impurity band to 
lower energies. We also notice that the impurity band in 
the ferromagnetic phase is broader than the one in the 
paramagnetic phase, confirming previous results [l2ll25|. 
This suggests that the lifetime of the "bound-state" holes 
is shorter in the ferromagnetic phase because additional 
scattering events that exchange holes between impurities 
are required to maintain magnetic order. Hence, we will 
refer to the impurity-band states as quasi-bound states. 

Next, we compare our density-of-states results with 
ARPES data in the paramagnetic phase of Gai-^Mn^As, 
x = 0.035 13] . For this doping Okabayashi et al. observe 
an impurity band already well separated but not very far 
apart from the main band. A rough estimate for J c of 
2.0 — 3.0 eV is the most suitable to describe the situa- 
tion observed in the experimental setup. In the following 



discussion we take J c = 3.0 eV, though this value is an 
overestimate of the exchange coupling, to further explore 
the physical consequences of our model and the behavior 
of the impurity band. 

Fig. |21 shows the decomposition of the DOS in terms 
of its J z components for J c = 3.0 eV. As expected, 
at T > T c all four components of J z contribute equally 
to the total DOS, the electronic system is unpolarized. 
At T = 0, however, we see that the impurity band is 
not fully polarized, as would be expected for the double 
exchange model. In addition to the dominant J z = +3/2 
component, components with J z = ±1/2 and —1/2 are 
also present. This is clearly a consequence of the strong 
spin-orbit coupling, which mixes the J z = ±3/2 state 
with J z = ±1/2 states. Previous DMFA studies[12( with 
coupling to only one carrier band captured the formation 
of the impurity band but were unable to address the effect 
of frustration on the carrier polarization. 

Now we explore the spectral function of our model for 
J c = 3.0 eV. Fig. 0(a) shows the spectral function in the 
paramagnetic phase. Figs.l^b) and (c) display the T = 
spectrum along the direction parallel and perpendicular 
to the average magnetization, respectively. Notice two 
main effects on the energy levels: the valence band quasi- 
particle states (shown in main panels) are renormalized 
and the impurity band (shown in the blown-up insets) ap- 
pears. Due to the localized nature of the impurity band 
states, their spectral weight extends over a large region in 
momentum space with typical values of the spectral func- 
tion reduced by two orders of magnitude in comparison 
with the quasiparticles peaks. 

The spectrum within the paramagnetic phase is 
isotropic, while it is obviously anisotropic in the ferro- 
magnetic phase. For J c = 3.0 eV all the quasi-particle 
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FIG. 3: Spectral functions near the center of the Brillouin 
zone for J c = 3.0 eV at (a) T > T c along any direction, 
(b) T — along the direction parallel and (c) perpendicular 
to the magnetization. Insets display a zoom of the impurity 
band region. Each spectral curve corresponds to a different 
value of the magnitude of the momentum k, as indicated at 
the right of the graphs, and the baseline of each spectrum is 
shifted proportionally to k for clarity. In the main panels the 
dashed curves represent the renormalized valence bands. For 
comparison, the valence bands of pure GaAs are also displayed 
(dotted curves). The scale in the insets is blown-up by a factor 
of 100 in (a) and 200 in (b) and (c). 



lines track the peaks, and quasi particles are well de- 
fined. As expected in the paramagnetic phase, the states 
corresponding to J • k = ±3/2 and ±1/2 remain degen- 
erate for all values of momentum. In this phase, the self- 
energy matrix is proportional to the identity [26|], pre- 
serving J • k as a good quantum number. The self energy 
just shifts the quasiparticle bands towards negative ener- 
gies due to the band repulsion between the quasiparticle 
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FIG. 4: Variation of the J z components of the impurity band 
spectral function for J c = 3.0 eV and T = ((a) J z = +3/2, 
(b) J z — +1/2, and (c) J z — —1/2) as the magnitude of the 
momentum k changes from (solid curve) to k — 0.25-7T along 
the direction parallel (dotted curve) and perpendicular (long- 
dashed curve) to the average magnetization. Notice that the 
three graphs have different scales. 



and the emerging impurity band. Also notice that the 
heavy and light quasiparticle bands are still degenerate 
at the r point. 

In the ferromagnetic phase, the quasi-particle lines 
split into four. Since the finite magnetization competes 
with the chiral nature of the holes, the quasi-particle lines 
no longer correspond to well-defined chiral states. This 
is confirmed by Figs. Ofb) and (c) where the curvature 
(related to the effective mass) of the quasi-particle bands 
depends on the direction of the momentum. For the di- 
rection parallel to the average magnetization the spec- 
trum of the quasiparticle band with J z = +1/2 is shifted 
towards the semiconducting gap by approximately OAeV 
while the others bands are pushed towards negative ener- 
gies. The positive shift in the J z — +1/2 band is mainly 
due to the Zeeman splitting. The band repulsion between 
this band and the impurity states is very small since the 
impurity band hardly includes J z — +1/2 quasiparti- 
cles with momentum parallel to the magnetization (see 
Fig. On the other hand the J z — +3/2 quasiparti- 
cle band shifts its spectrum by w ~0.2eV due to strong 
repulsion with the J z = +3/2 states at the impurity 
band. However, for the direction perpendicular to the 
average magnetization the quasiparticle peak centered at 
u> w OAeV has J • k = +3/2 character. Since the impu- 
rity band on the perpendicular plane is formed mainly 
by J • k = ±1/2 states, the J • k = +3/2 quasiparticles 
do not suffer band repulsion. 

Now examine the impurity band spectra in the insets. 
In the paramagnetic phase the impurity band does not 
show significant variation with k, indicating strong lo- 
calization of the bound-state holes. In the ferromagnetic 
phase, the variation with k, although small, confirms that 
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the bound-state holes are more mobile. Typical fillings 
leading to the highest critical temperatures correspond 
to values of the chemical potential inside the impurity 
band. Therefore, the transport properties rely on the 
impurity band rather than on the quasi-particle bands. 

As demonstrated by Fig.^for the ferromagnetic phase, 
the impurity-band spectral function along the direction 
parallel to the magnetization is predominantly composed 
of J z = +3/2 holes, whereas along the direction per- 
pendicular to the magnetization, it is a mixture of J z — 
+3/2,-1-1/2, and —1/2 states. This result suggests that 
the carrier polarization may be optimized by driving the 
current along the direction parallel to the magnetization. 
However, this may come with a price, since the carrier 
mobility is lower for the J z = +3/2 heavy holes. Due to 
the mixing of heavy and light holes, the direction per- 
pendicular to the magnetization may have higher mobil- 
ity but lower carrier-spin polarization. Also notice that 
the J z = —1/2 states participating in the local screening 
at the impurity band mostly display finite perpendicular 
momentum. This is the configuration most energetically 
favorable to avoid an exchange penalty while fulfilling 
spin-orbit constraints. 



In conclusion, we have calculated the spectra of the 
renormalized valence bands and the impurity band of 
Gai-^Mn^As within the DMFA. We compare our results 
with existing ARPES data for the paramagnetic phase 
|13| . From the anisotropy of the impurity band in the 
ferromagnetic phase, we predict that the direction paral- 
lel to the magnetization will produce the most polarized 
spin current, whereas the perpendicular direction may 
display higher conductivity with lower polarization. It 
would be interesting to be able to compare our results in 
the ferromagnetic phase with additional ARPES data. 
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